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1 Introduction 


The fragmentation of quarks into hadrons is still not completely understood. In fact, 
the final stages of hadronization involve low scale nonperturbative processes in QCD in 
Minkowski space together with hadronic wave functions which are not yet calculable. The 
nucleus may help to understand the space time evolution of a parton as the nucleons inside 
the nucleus play the role of very nearby detectors of the propagating object. 

Deep inelastic lepton-nucleus scattering has the advantage that the lepton transfers a 
well-defined energy v to the struck quark propagating through cold nuclear matter. Due 
to factorization in deep inelastic scattering (DIS), the semi-inclusive cross section can be 
described by the product of a parton distribution function (PDF) with a fragmentation 
function (FF) and the hard scattering cross section. Fig. 1 shows a schematic diagram 
of semi-inclusive deep inelastic lepton scattering on a target, and the definitions of the 
four momenta of the particles involved in the process. In the naive parton model, the 
probability qj(x) that a quark of flavor / with momentum fraction x is present in the 
target is multiplied with the probability D^(z) that it hadronizes into a definite hadron 
h which carries a momentum fraction z of the quark. 

While the inclusive deep inelastic scattering of high energy leptons on nuclei has been 
used to measure the quark distributions q/(x) in nuclei, semi-inclusive DIS (SIDIS) can 
be used to study the medium modifications of quark fragmentation functions Dj(z). In 
SIDIS, besides the scattered lepton l', the leading hadron h formed from the struck 
quark is detected with energy Eh = zv in the target rest frame (see Table 1 for a 
list of kinematic variables). Experimentally, hadron multiplicities are obtained from 
normalizing the SIDIS yield N h to the DIS yield N DIS \ (1 /N DIS )(dN h (x, z)/dz) ~ 
(52/ e 2 f q f (x)D h f (z)/52 f e /?/( a; ))) where the sum over all the quark of a flavor / and charge 
Cf is performed. It has been shown that, after integrating over a broad range of x, the 
multiplicities are a good approximation to the fragmentation functions (see Ref. [1] and 
references therein). Thus, a measurement of the hadron multiplicity in nuclei is a good 
tool to measure possible medium modifications of quark fragmentation functions. 

The experimental results on semi-inclusive leptoproduction of hadrons from nuclei 
[2,3] are usually presented in terms of multiplicity ratios between nuclear (A) and 
deuteron (D) targets, as functions of z and v respectively: 


1 dN\ / 1 dN h D 
~ N£ is dz / N£ IS dz 
1 dN\ / 1 dN h D 
~ N£ is dv / n£ is dv 


( 1 . 1 ) 

( 1 . 2 ) 


In absence of nuclear effects, the ratio R h M would be equal to 1, and experimental 
results show that this is the case at high transfer energy v [2], 

The HERMES experiment [3] has brought new insight into the study of the hadroniza¬ 
tion process, since it has access to lower energy transfers z/, where the hadronization is 
supposed to occur inside the nucleus. It has also extended the range of fractional momenta 
to larger values of z. In principle, a comparison with the EMC results [2], which involve 
much higher transfer energy, would allow us to separate the physics of gluon radiation 
from the absorption process of the produced hadron. 
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Four-momenta of the particles 


Particle 

In lab. 

frame 

Incident lepton 

k = 

C E,k ) 

Scattered lepton 

k! = 

(E',k’) 

Target nucleon/nucleus 

P = 

(M,0) 

Virtual photon 

q = 

(y,v) 

Produced hadron 

p = 

( E h ,p) 


Figure 1: Semi-inclusive hadron production in deep inelastic scattering on a target T in the pQCD 
factorization approach. Parton distribution functions (PDF) and fragmentation functions (FF) represent 
the non-perturbative input. Four-momenta of the particles involved in the process in the laboratory 
frame are defined in the table on the right. 


There has been an intensive discussion on energy loss of partons in a cold and hot 
nuclear environment [5], which then leads to modified fragmentation functions in nuclei 
or in deconfined quark matter [6]. Recent data on high transverse momentum hadron 
produced in Au-Au collisions at RHIC have truly opened up the domain of hard scattering 
and jet-quenching in heavy-ion collisions [7]. We would like to address this problem from a 
different perspective. Already in the vacuum, the struck quark loses substantial energy in 
the gluon cascade; the question arises how this cascade process is modified in the nuclear 
medium. 

The goal of this paper is to describe the R\j measurements [2,3] 1 in terms of modified 
fragmentation functions supplemented by nuclear absorption. We describe the modifica¬ 
tion of the structure function and fragmentation function in the framework of rescaling 
models [8,9,11]. These models have been successfully applied to the interpretation of 
nuclear structure functions. In this paper, we extend the rescaling model also to the 
medium modification of the quark fragmentation functions. For the nuclear absorption, 
we work out two specific fragmentation models based on work by Bialas-Gyulassy [12] 
and Bialas-Chmaj [13]. Applications of the above mentioned nuclear absorption models 
to DIS on nuclear targets may be found in Refs. [14] and [15]. 

Recently, other “QCD inspired” analysis such as the gluon-bremsstrahlung calculation 
for leading hadron production [16], the higher-twist pQCD computation of Ref. [17] and 
the effect of parton energy loss of Ref. [18] have been applied to hadron production in DIS 
on nuclei. Computations for meson production in DIS on nuclei in terms of stationary 
string model are also presented in Ref. [19]. 

The outline of the paper is as follows: in Sec. 2, we describe the modification of 
the fragmentation function by gluon radiation of the struck quark; in Sec. 3, we discuss 
the formation length and in Sec. 4, the nuclear absorption of the produced hadron. The 
comparison between the model predictions and the experimental data for charged hadrons 
is made in Sec. 5. The flavor dependence of the hadronization process is calculated and 
compared with data in Sec. 6. Conclusions are reported in Sec. 7. 

1 The SLAC data [4] will not be considered in this paper since the semi-inclusive distributions have 
been normalized to the target density rather than to the rate of inclusive events. Therefore they can not 
be directly compared to EMC and HERMES data. 
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Kinematic variables 


Variable 

Covariant 

Lab. frame 


Q 2 = 

—q 2 ~ 

2Mxp 

Negative four-momentum squared 
of the virtual photon. 

V = 

q-p — 

Vp ^ 

E' - E 

Energy of the virtual photon in lab. 
frame. 

X = 

-q 2 

2 Pq 

Q 2 

2 Mu 

Bjorken scaling variable (fraction of 
P carried by the struck parton in 
the Breit frame). 

2 = 

•a re 

‘'altl 

II 

Eh 

V 

Fraction of the virtual photon en¬ 
ergy carried by the hadron. 

y = 

ariQ 

ii 

V 

E 

Fraction of the incident lepton en¬ 
ergy transferred to the target. 

w 2 = 

(p + q) 2 = 

M 2 + 2 Mu - Q 2 

Invariant mass squared of the total 
hadronic final state. 


Table 1: Definitions of the kinematic variables used in the paper. 


2 Modification of parton distribution and fragmen¬ 
tation functions by gluon radiation. 

It has been shown (for a recent review see Ref. [20]), that a simple folding of nucleon 
structure functions with the momentum distribution of nucleons in the nucleus is insuffi¬ 
cient to reproduce the nuclear structure functions. Nuclear structure and fragmentation 
functions certainly depend on the density of baryonic matter and on the amount of mo¬ 
mentum sharing induced by the close packing of nucleons. A good description of nuclear 
structure functions has been achieved in the so-called deconfinement models, based on 
the hypothesis that quarks in bound nucleons have access to a larger region in space than 
in free nucleons [8-10]. 

Both parton distribution functions and fragmentation functions depend on the virtu¬ 
ality Q 2 of the DIS process. Their adjustment to the physical scale Q 2 takes into account 
all radiated gluons before and after the photon-quark interaction in the leading logarithm 
approximation. The long wavelength spectrum of gluons extends farther into the infrared 
toward low Q 2 oc 1/A 2 where A is the confinement scale. Deconfinement models assume 
a larger confinement scale A^ in nuclei, compared with the confinement scale Ao in free 
nucleons: 


An > Ao • 

Therefore, in nuclei gluon radiation would be affected by the assumed deconfinement of 
color. 

Take, e.g., a parton distribution function and consider a quark which carries a mo¬ 
mentum Qq when it is confined on a scale Ao- If the scale changes to An it carries a 
corresponding momentum Qa■ If the free and bound nucleon would be characterized by 
these scales alone, the structure function of the bound nucleon would be related to the 
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structure function of the free nucleon by replacing the measured Q 2 by (Aa/Ao) 2 Q 2 in 
the arguments of the structure and fragmentation functions. Taking into account also 
the running of the QCD coupling, i.e. the existence of the QCD scale A qcd, the correct 
rescaling factor is [9,10]: 

UQ 2 ) = (^) ^ , (2.i) 

where a s (Q 2 ) is computed at leading order with Aqcd = 200 MeV and four quark flavours. 
The real scale factor £a(Q 2 ) considers that the low lying modes couple with a coupling 
constant 6t s . Effectively, the DGLAP evolution of the nuclear structure function covers a 
larger interval in momentum compared with the corresponding functions in the nucleon 
at the same scale Q, and the nuclear modified distribution function qj- for a quark of 
flavour / reads: 


qj(x,Q 2 ) = qf(x,U(Q 2 )Q 2 ) , 


( 2 . 2 ) 


where qj is the corresponding distribution function in a free nucleon. For consistency 
it is necessary that partial deconfinement in nuclei does not only modify the parton 
distribution functions but also the fragmentation functions [11]. A procedure similar to 
the above one for the fragmentation function gives 

D>} A (z,Q 2 ) = D h f (z,U(Q 2 )Q 2 ) , (2.3) 

where D (5 is the fragmentation function of a quark of flavour / into a hadron h and !/'} ' 
is the nuclear modified fragmentation function. Of course such a simple adjustment of 
evolution scale is only a simple prescription to mimic the effects of the nucleus on the 
light-cone wave function of the bound nucleon. The rescaling relations (2.2) and (2.3) are 
expected approximately to be valid for 0.1 < x < 0.6 and z > 0.1. 

We use two models for the computation of the scale factor. The maximal deconfine- 
ment model (MD) [8,10] assumes the onset of “colour conductivity” in nuclei and takes 
Aa = Ra, where Ra is the nuclear radius, so that 


Aa _ Ra 
Aq Rp 


(2.4) 


Since in this model the low lying modes can become very soft a s is frozen at a s = 0.54. 
The partial deconfinement model (PD) [9] assumes the deconfinement scale Aa to be 
proportional to the overlap of nucleons inside the given nucleus, and a s = a s (fi\), where 
Ha = j^H o and Ho = 0.66 GeV 2 is the starting scale for DGLAP evolution in the free 
nucleon. A recent analysis [21], based on the experimental results [22] on the nucleon 
form factor in nuclei, found that the amount of deconfinement predicted by Ref. [9] is too 
large. In particular Ref. [21] gives an upper limit for (Aa — A 0 )/A 0 ~10% on iron, while 
Ref. [9] predicts a 15.3% enhancement. Therefore following Ref. [21] we reduce the Aa 
computed in Ref. [9] by a universal factor: 


(\A 

Wo 



0.10 

0.153 



CJRR 


(2.5) 
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2 Dr 14 N 7 20 Neio 56 Fe 26 63 Cu 29 84 Kr 36 131 Xe 54 

PD model 

A a Mo 

VU(Q = 1.5 GeV) 

1.010 1.071 1.068 1.100 1.101 1.110 1.118 

1.014 1.110 1.106 1.164 1.165 1.182 1.195 

MD model 

A a Mo 

VCAiQ = 1-5 GeV) 

2.46 2.99 3.53 4.42 4.62 4.88 5.68 

3.65 4.82 6.13 8.47 9.02 9.76 12.3 


Table 2: Values of the deconfinement ratio A^/Ao and scale factors used in partial deconfinement 
(PD) model (see Eq. (2.5)) and maximal deconfinement (MD) model (see Eq. (2.4). 


In Table 2 we show the values of the scale factor M predicted by the aforementioned 
maximal deconfinement (MD) and partial deconfinement (PD) models for several nuclei. 

The effect of the rescaling on the fragmentation function D* + (z, C,aQ 2 ) [23] for different 
nuclei is shown in Fig. 2 both for the partial and maximal deconfinement models. Note 
that the MD model assumes a much larger scale factor, which results in a larger nuclear 
modification of hadron production especially at high z. One also sees that in the partial 
deconfinement model the effect of rescaling on the deuteron is much smaller than on heavy 
nuclei. This is different from the maximal deconfinement model, where only the nuclear 
radius matters and rescaling effects gradually vary between light and heavy nuclei. It is 
also worth to point out that the fragmentation function ratios shown in Fig. 2 increase in 
the low -2 region up to values larger than unity approximately at 2 ~ 0.2 for all nuclei in 
both rescaling models. This increase, common to all fragmentation functions with rescaled 
Q 2 , comes from enhanced gluon radiation at small 2 contributing to hadron production. 




Figure 2: Ratio of the rescaled fragmentation functions u — > 7r + to the standard ones at leading 
order [23] for Q 2 = 2.25 GeV 2 as a function of z, in the case of partial deconfinement rescaling (PD) 
and maximal deconfinement rescaling (MD). 

We calculate the multiplicity ratio of Eqs. (1.1) and (1.2) by using the rescaled PDF 
(2.2) and the rescaled FF (2.3) in the leading order pQCD computation of N ^ IS , dN\/dz 
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and dN\jdv\ 


1 dN\ 
N° IS dz 

1 dN\ 
N® IS dv 

a iA 


a 


<7 


1 

7a 

1 

1a 


/ dx dv E ej q f {x,UQ 2 

*exp. cuts j 


/ dxdz E e 2 f Qf(x,U Q ' 2 

*exp. cuts ^ 


/ dx dv E e )q f (x,UQ 2 

' exp. cuts j 


da lq 

dxdv 


t dv D '^) 

(2,6) 

L, d ' ( ^ q2) 

(2,7) 

q 

(2,8) 


Here da £q /dxdv is the differential cross-section for lepton-quark scattering computed at 
leading order (LO) in perturbation theory [24]: 


do** 

dxdv 


Mx 


4na 2 (Q 2 ) 

Q 4 


[i+(i-i/) 2 ]. 


Note that the experimental acceptance is explicitly accounted for in the integration limits 
in Eqs. (2.6)-(2.8). The experimental acceptances of the EMC and HERMES experiments 
are given in App. A. Nuclear isospin asymmetry is taken into account by using the 
averaged quark distribution function 




The quark distribution function in the proton and in the neutron, respectively, and qj, 
are related by isospin symmetry. In the numerical computations we use the leading order 
GRV98 parton distribution functions [25]. Since the PDF’s appear both in the numerator 
and in the denominator of Eqs. (2.6)-(2.7), the theoretical uncertainty resulting from a 
different parametrization of the distribution functions is minimal. For the same reason, 
the effect of shadowing corrections to PDF’s in the small-x region is negligible (the lowest 
x for both EMC and HERMES experiments is x ~ 0.02). We checked that by using the 
EKS98 parametrization [26] of PDF’s with shadowing corrections, the multiplicity ratios 
(1.1) and (1.2) change by less than 1%. 

Furthermore, we analyzed charged hadron production with both the Khnicl-Kramer- 
Pbtter parametrization of the FF at leading order (KKPLO)[27] and the leading order 
Kretzer’s parametrization (KLO) [23]. The different parametrization of the FF produces 
a 1% difference in the multiplicity ratios. In the rest of the paper we shall use the 
KLO parametrization of the FF since it gives charge-separated fragmentation functions 
for pions and kaons, which are necessary for a comparison of the theoretical model with 
the HERMES data. Note, however, that by isospin symmetry arguments it is possible 
to relate charge-averaged FF’s, like KKPLO, to charge-separated ones at the cost of an 
additional parameter [28]. 

Due to Lorentz dilatation, at large values of v the hadrons are expected to form 
mainly outside the nucleus, so that the effect of reinteractions of the hadron with the 
nucleus arc minimal. For this reason, the EMC data on hadron production in high energy 
DIS scattering on copper are ideal for a comparison of the partial and the maximal 
deconfinement model. The results for the z- and //-dependence of the multiplicity ratio 
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Figure 3: Theoretical multiplicity ratios Rm as functions of z and v in the PD (solid lines) and 
MD (dashed lines) models are compared with EMC data [2] on Cu. The theoretical curves between 
10 GeV < v < 30 GeV are computed with the experimental acceptance corresponding to a beam 
energy of 100 GeV (see App. A). 


are shown in Fig. 3 and compared with EMC data [2], The maximal deconfinement model 
(dashed lines) underestimates the data over nearly the whole range in z and v. On the 
contrary the partial deconfinement model (solid lines) better reproduces both the z and 
v experimental distributions but overestimates the data at low v and z, thus suggesting 
that additional processes contribute in this region. As the maximal deconfinement model 
overpredicts the effect of deconfinement in the nucleus, in the following we shall consider 
only the partial deconfinement model. 

3 Formation length distribution in the Lund model 

The results presented in Fig. 3 show that the partial deconfinement rescaling model overes¬ 
timates the multiplicity data slightly at small energy transfers v and at small z. Therefore, 
we have to consider in addition the formation of the hadron and its subsequent interaction 
in the nuclear medium. As the hadron formation length decreases for lower v, the effect 
of nuclear interaction becomes more relevant in the kinematic region of the HERMES 
experiment, and the effect is amplified in a heavy target as the formation length is com¬ 
parable with the nuclear radius. For a realistic calculation, a good understanding of the 
formation mechanism is necessary. 

In this paper we follow the LUND model for the fragmentation process [29, 30] and we 
apply the hadron formation and interaction discussed in Ref. [12] to the DIS process. A 
picture of the space-time development of the fragmentation process in the Lund model is 
presented in Fig. 4. The process begins at t=0 and y=0 when the quark q is ejected from a 
nucleon by the virtual photon . The quark propagates in the positive longitudinal direction 
y and a colour string is formed between the quark and the target remnant. The string 
maps out the lightly shaded space time area. It breaks at points C t into smaller pieces 
due to quark-antiquark pair creation. When a quark moving in the positive direction 
meets an antiquark moving in the negative direction, they interact and form a so called 
“yo-yo” state, which is identified with a final state meson. This happens at points P % in 










Figure 4: Formation lengths of the pre¬ 
hadron and hadron are shown in the Lund 
string model. The shaded area represents the 
world-sheet swept by the string. Produced 
hadrons, hi, are ordered according to their 
rank i (here i = 1,2,3 and 3 hadrons are pro¬ 
duced). At points Ci prehadrons are created 
and at points Pi the hadrons hi are formed. 
The formation lengths Z* and lh at the bottom 
of the figure refer to the second-rank hadron 
h 2 . 


the figure, where the created hadron is denoted by hi. We shall denote by Zi the fraction 
of the energy it carries. Hadrons are ordered according to their rank i [29]. Note that the 
first-rank hadron is always created at the end of the string, whose length is in our case 



where Ka is the string tension in the nuclear environment. If partial deconfinement occurs, 
also the string tension is rescaled in the nucleus: 

= i ( 3 - 1 ) 

where k ~ 1 GeV/fm is the string tension in the vacuum. Since the string tension is the 
physical quantity that sets the confinement scale, a larger confinement scale corresponds 
to a smaller string tension. 

Let us focus our attention on the h 2 hadron in Fig. 4. There are two relevant lengths 
for the fragmentation process: 

(%) prehadron (or constituent) formation length /* at which the first constituent of h 2 
is created. 

(ii) the hadron (or yo-yo) formation length G < L at which the hadron h 2 is formed. 

In the Lund model these two lengths are related in the following way [12, 30]: 

lh — l* + ^ 2 L , (3-2) 

At fixed z they both increase linearly with the virtual photon energy v. However, as 
functions of z they behave rather differently, especially at z—*l, where /*—>0 and lh~*L 
(see Fig. 5). Therefore, we have two distinct scenarios of fragmentation, whether the 
prehadronic state interacts with the nucleus, or not. Analysis of the data on hadron- 
nucleus and lepton-nucleus collisions strongly suggest that the prehadron behaves as an 
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Figure 5: The average normalized pre¬ 
hadron formation length (Z*)/L computed in 
standard Lund model Eq. (3.4) with C=0.3 
(solid line), compared with the one calcu¬ 
lated from Eq. (3.5) used in the EMC analy¬ 
sis (dotted line). The normalized hadron (or 
yo-yo) formation length ( lh)/L is shown as a 
dashed curve. 


hadron-like object long before final state hadrons actually appear [12,31]. Also an analysis 
of EMC and SLAC data confirms the importance of the prehadron interactions [15]. 

In App. B we derive the probability distribution V * that a prehadronic state is formed 
at a formation length y with momentum fraction z. In the case that the prehadron is 
composed of q = u,d,s flavours, the distribution reads 


'P*(y, z, L) = 


zL 


y 


y-zL [(y + zL)(l -z) 
x < S[y - (1 


c 


1 - 4 - C 

z)L\ + -±-0[( 1 


y 


zL 


z) L ~y\ \ %], 


(3.3) 


where C ~ 0.3 is the only parameter entering this expression. We can then compute the 
hadron average formation length (If), which we identify with the prehadronic formation 
length, (If) = (U)(z,L) = f dyyV*(y, z, L): 


(If) - 




1+C 1 


(1 — z) z L , 


(3.4) 


where 2 F 1 is the Gauss’ hypergeometric function [32], Note that as z—*l the average 
formation length (If )—->(1 — z)- behaves similarly to the formation length suggested by 
the gluon bremsstrahlung model of Ref. [16]. Indeed, quantum mechanics yields a large 
energy for an intermediate state consisting of a quark and a gluon with z and 1 — z 
momentum fractions when z—>l. Consequently, if the hadron takes almost all the energy 
(z— >1) then the quark-gluon system is short lived and must hadronize almost immediately 
[33]. 

In Fig. 5 we show the average formation length computed with Eq. (3.4) and C = 0.3 
as in the standard Lund model (solid line). This is the formula we use in the next 
sections. For comparison (dotted line) we show also the average formation length used in 
the analysis of the EMC data. [34], 


(If) 


Tn(l /z 2 ) -1 + z 2 
1 — z 2 


zL , 


(3.5) 
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Figure 6 : Nuclear attenuation of a hadron 
ft: the virtual photon 7 * interacts with a 
quark q at a longitudinal coordinate y\ the 
quark turns into a “prehadronic” state ft* at 
position y' and the hadron ft. is formed at 
y". Each state interacts with the surround¬ 
ing nucleons with a cross-section a q , cr* and 
ah, respectively. 


which overestimates up to 30% the Lund model result. The above formation length is 
obtained in the Lund model with a non-standard choice of parameters, as described in 
App. B. Finally, the contribution of the hadron (or yo-yo) formation lenght % hr the 
standard Lund model is shown with a dashed curve. 


4 Absorption of the produced hadron in the nucleus 

The spatial evolution of a quark q created by the virtual photon 7 * at a longitudinal 
position y is shown in Fig. 6 in all intermediate stages. In the first stage the quark 
propagates to the position y' > y where a prehadronic state h* is formed. In the second 
stage the final hadron h is created at the point y" > y'. 

The quark, the prehadronic state and the final hadron propagate through nuclear 
matter and interact with the surrounding nucleons. As a result, each of the three states 
may undergo inelastic interactions and/or lose longitudinal energy through elastic rescat¬ 
terings. In a first approximation we may assume that the final hadron with a value of 
z corresponding to the fragmentation process in vacuum will be observed provided that 
none of the three propagating states has interacted with the nucleus. By defining the 
nuclear absorption factor A/q (z, u) as the probability that neither the q, h„ , nor h have 
interacted with a nucleon, the multiplicity ratios of Eqs. (2.6) and (2.7) are modified as 
follows: 

1 HN h 1 f _ drr eq 

n dis dz A = Y, e / q f( x ’ ^ Q 2 ) /TAfi) D f^ ^ g2) Na{z ' ’ t 4 ' 1 ) 

1 HN h 1 ft 

^ J^^)v^U Q2 )^D)(z,UQ 2 )N A (z,u) . (4.2) 

While the fragmentation functions Dj are sensitive to the virtuality Q 2 of the fragmen¬ 
tation process in the medium, the nuclear absorption factor A (4, which depends on z and 
v, is sensitive to the hadron energy in the rest frame of the nucleus. 

We will investigate two models for the computation of the nuclear absorption factor: 
the Bialas-Gyulassy (BG) model, [12], and the Bialas-Chmaj (BC) model, [13]. These 
models require the probability distribution function of the hadron formation length and 
the average formation length discussed and derived in Sec. 3. In order to reduce the 
number of the parameters of these models we set the quark-nucleon cross-section cr q = 0, 
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since the analysis of the DIS data on R\j reported in [ 2 , 3, 34] does not indicate a significant 
interaction of the quark with the nuclear medium. Moreover, we do not give the formulae 
for a separate evolution of the prehadron to the full hadronic size in the main text. 
Then, the prehadronic state h* and the hadron h are treated effectively as a single object 
h* = h created at y' which propagates in the nucleus with a single cross-section cr* = cr h . 
The effect of distinguishing the prehadronic state by a different cross section <v* ^ 07 , is 
discussed in Sec. 5 and 6 , and the formulae needed in the computations of the absorption 
effect are discussed in App.C. 

4.1 Bialas-Gyulassy (BG) model 

In the Bialas-Gyulassy (BG) model [12] the nuclear absorption for a one-step fragmen¬ 
tation process with cr* = a h depends on three parts: the nuclear density profile p A (b,y) 
normalized to unity as a function of impact parameter b and longitudinal coordinate y\ 
the probability distribution V* (y ' — y, z, L) that the prehadron is formed after the length 
y' — y; and the probability [Wo^O ]" 4-1 that the hadron created at y' does not interact 
with the nucleus. The nuclear absorption factor has the following form: 

/ roo /»oo 

d 2 b dyp A (b,y) dy'V*(y' - y, z, L) [W 0 (y')] A_1 » (4.3) 

J —00 J y 

where 

/»oo 

Wo(y > ) — 1 — a h / dy*p A (b,y*) . (4.4) 

Jy' 

For 7A we will use the probability distribution (3.3) discussed in Sec. 3. 

4.2 Bialas-Chmaj (BC) model 

The Bialas-Chmaj (BC) model [13] is a simpler version of the BG-model for including 
nuclear absorption. In the BC model the nuclear absorption factor depends on the prob¬ 
ability P h {y',y ) that the hadron h* = h is formed at a distance y' — y from the 7*<7 
interaction point, 

P h (y'-y) = i-e-W-vVVr) , (4.5) 

and on the survival probability S A {b,y) A ~ l of this hadron, where for a one step fragmen¬ 
tation process with cr* = ah one has 

/»oo 

S A (b,y) = l-a h dy' P h {y' - y)p A {b,y > ) . (4.6) 

Jy 

The nuclear absorption factor is then written as: 

/ /»oo 

d 2 b / dyp A (b,y)[S A (b,y)] A 1 , 

J —OO 

where the average formation length (l F ) = {l F ){z, v) we will use has been derived in Sec. 3, 
Eq. (3.4). 
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Figure 7: Isospin averaged hadron- 
nucleon total cross section versus the 
hadron energy Eh = Ei a b , in the case 
of krypton nucleus, for positive and 
negative pions and kaons, and protons 
and antiprotons. The charged hadron 
cross-section ah has been computed as 
weighted sum of individual cross sec¬ 
tions for pions, kaons, protons and an- 
tiprotons. 


To understand the BC model and its relationship to the BG model, Eq. (4.3), we can 
approximate p^(p, b) with a uniform nuclear density Pa(v, b) = Po/A inside a hard sphere. 
Here p 0 is the nuclear matter density. In this approximation the integrations in Eq. (4.6) 
may be carried out explicitly and yield 


A/f c) 


1 _ {R( b) - (j, _ ((f)) 0 (R(b) -y)e(R-b)} 


(4.7) 


Here R{b) 2 = V R 2 — b 2 and the effective formation length (l F ) is defined as: 



(If) 


e -(R(b)-y)/(l F ) 1 


(4.8) 


It is then easy to see that Eq. (4.7) can be derived from the Bialas-Gyulassy formula (4.3) 
if we choose 

V,(y'-y)=S[(y' -y)-{lf)} . 

In this sense, the BC model for nuclear absorption is an effective model for absorption 
which describes the fragmentation process by modifying the average formation length ac¬ 
cording to Eq. (4.8), and letting the hadron to be formed suddenly after a well dehned 
formation length, instead of assuming a probability distribution in formation lengths. 
Nonetheless, as we shall discuss in Sec. 5 and Sec. 6, the BC model gives a useful phe¬ 
nomenological description of nuclear absorption processes. 


5 Numerical results for multiplicity ratios in nuclei 

In this section we discuss the results of our model for semi-inclusive charged hadron 
production in nuclei including effects of perturbative gluon radiation and absorption, and 
compare them to EMC and HERMES experimental data. 

In the computation of the nuclear absorption factor, we use as a density for the 
deuteron the sum of the Reid’s soft-core S- and D-wave functions squared [35]. For 
heavier nuclei we use a Woods-Saxon density with radius Ra = 1.12H 1 / 3 — 0.86H -1 / 3 fm. 
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As discussed in Sec. 4, the fragmentation of the quark is considered without differenti¬ 
ation between the prehadron-nucleon and hadron-nucleon cross sections. For the hadron- 
nucleon cross sections, we take the isospin averaged hadron-proton Uhp and hadron- 
neutron < 7 hn total cross sections: 

a h (E h ) = [. Za hp (E h ) + (A - Z)a hn (E h )]. (5.1) 

As the total hadron-nucleon cross sections include also the elastic ones, the correspond¬ 
ing results represent an upper limit with respect to the nuclear absorption effect in the 
multiplicity ratios. The energy dependence of (Thp and ah n [36] has been explicitly taken 
into account, and is at variance with respect to the original models [12,13]. The en¬ 
ergy behavior of the total hadron-nucleon cross section is shown in Fig. 7, in the case 
of a krypton target, for positive/negative pions and kaons. The corresponding charged 
hadron cross section has been derived as weighted sum of individual cross sections for 
pions, kaons, protons and antiprotons. The weighted sum has been performed by taking 
into account the measured relative yield of various hadrons reported by the HERMES 
experiment in the 3-23 GeV region [3] as this region is very sensitive to the variation of 
the hadron-nucleon cross section. 

As the model predictions are sensitive to the average kinematic variables used in the 
computations, in Fig. 8 we compare our LO theoretical computations of the average is, 
z and Q 2 with the corresponding experimental values. The model reproduces reasonably 
well the experimental (Q 2 ), ( z ) and (u) behaviour. The step structures found in the LO 
computation of (is)(z) and (. z)(is ) are due to the experimental kinematic cuts shown in 
Fig. 12 of App. A. The differences compared with the experimentally observed smoother 
dependence of (is) and (z) may come from NLO corrections or higher-twist terms. There¬ 
fore, in the computation of the multiplicity ratios we replaced the arguments of Eq. (4.1) 
by experimental values when available: 

Qf{x,U * Q.fip'i (Q")exp(z)) 

d ^. uq 2 ) —* o}(u<(e’U(!)) 

Na{z,v) -> (v)e X p(z)) , 

and in Eq. (4.2) 

qf(x,UQ 2 ) —♦ qf(x,U(Q 2 )e X P H) 

D h f (z,UQ 2 ) — D h f (z,U(Q 2 )e, M) 

J\f A (z,is) -» J\f A ((z) exp (is),is) . 

The effect of using experimental averages instead of the LO-computed ones by Eqs. (4.1) 
and Eq. (4.2), maximally yields 2-3% variation in the presented multiplicity ratios. The 
multiplicity ratios for positive and negative charged hadrons are calculated by using the 
fragmentation functions for hA in Kretzer’s parameterization, with the meson formation 
time in the Lund model, Eqs. (3.3) and (3.4), and with the average cross-section shown in 
Fig. 7. Since charged hadrons are dominated by meson production, this is a good approx¬ 
imation to the sum of individual hadron fragmentation yields multiplied by individual 
absorption factors. In Sec. 6 we will also calculate multiplicity ratios for pions and kaons 
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Figure 8 : Average theoretical values of kinematic variables computed with the LO formulae (2.6) and 
(2.7) are compared with the experimental kinematic variables at HERMES for 7r and K production on 
Kr [37]. Left: {v) and ( Q 2 ) as a function of 2 . Right: ( z ) and ( Q 2 ) as a function of v. 


separately. We leave out an extra calculation of baryon and antibaryon production, since 
such a calculation would require to know the contribution of target fragmentation and 
to extract separate fragmentation functions for p and p from the charge averaged KKP 
parameterization. Furthermore a different formation mechanism is possible for baryons 
or antibaryons, see Eqs. (B.15) and (B.16). 

The theoretical calculations of the multiplicity ratio for charged hadrons are presented 
in Fig. 9 and compared with the EMC and HERMES data for different nuclei. In the 
figure, the three different theoretical curves correspond to calculations with rescaled frag¬ 
mentation functions alone (dashed curve), with rescaling plus nuclear absorption in the 
BC model (solid curve), and with rescaling plus absorption in the BG model (dotted 
curve). 

In the case of BG absorption with the probability distribution of Eq. (3.3), the general 
shape of the curve deviates from the experimental behavior and the calculation under¬ 
estimates the multiplicity ratio both at small and large z. In the case of BC absorption 
there is a better agreement with the experimental behaviour. Note that in this case we 
used the average formation length of Eq. (3.4) where we fixed the effective string tension 
k-bc — 0.4 GeV/fm in order to agree with the experimental data. As already observed in 
Refs. [2,34] and [38] the Bialas-Chmaj model with k — 1 GeV/fm overestimates by far 
the absorption effects, while, with a reduced effective string tension, a reasonably good 
description of the data is obtained inside the statistical and systematic errors. 

The observed differences between the data and the model predictions in the low- 2 ; 
region can be interpreted in several ways. In the described models for hadron absorption, 
a hadron which interacted with a nucleon is simply removed without accounting for a 
possible lowering of 2 clue to the interaction. A modification of the string fragmentation 
mechanism to take this into account was proposed in Ref. [14]. Moreover, the differences 
at low 2 can be ascribed to target fragmentation protons, which can give a significant 
contribution to the multiplicity in this region. Also the effect of different prehadron and 
hadron cross sections <r* ^ oy may play a role. This has been investigated by fixing 
the prehadronic state cross-section cr* = 0.5 ov The corresponding results, presented in 
the bottom right plot of Fig. 9, show an enhanced multiplicity ratio, but the shape of 
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Figure 9: Theoretical multiplicity ratios of charged hadrons for Cu (EMC), Kr and N (HERMES) 
targets as function of 2 compared with the data [2,3]. Dashed lines show the predictions without 
absorption and only rescaling, according to Eqs. (2.6)-Eqs. (2.8); solid lines give the calculation 
with rescaling and absorption according to Eqs. (4.6) and (3.4) for the BC model; dotted lines 
with rescaling plus absorption according to Eq. (4.3) and (3.3) for the BG model. In the bottom 
right plot, the solid and dotted curves represent the result of the computations with er* = 0.5 07 *. 
The error bars represent the statistical uncertainty only. Systematic uncertainties for the HERMES 
nitrogen and krypton data are 3% and 3.5% respectively [3,37]. 


the BG predictions still disagrees with the data and the BC predictions overestimate the 
experimental results. Therefore, we will use the BC model with prehadronic cross section 
equal to the hadronic cross section to predict R\j in the region 0.2 < z < 0.9 also for 
different fragmentation products and for different targets (see Sec. 6). 

As we can see in Fig. 9, the effect of the absorption is less important at EMC in Cu than 
the absorption effect at HERMES in Kr, which is a nucleus of similar size, due to the larger 
energy transfer v in the EMC data ((z/)=62 GeV) compared to HERMES ((z/)~ 12 GeV). 
In the HERMES kinematic range both rescaling and absorption contribute. Absorption 
becomes the dominant effect in Kr tending to mask the rescaling effect. In the case of 
the N nucleus the theoretical result underestimates the experimental data. This may be 
due to an overestimate of the rescaling effect in light nuclei, which have a large surface to 
volume ratio. 

It is worth to point out that the accuracy of the HERMES data may distinguish 
different models for the space-time development of the fragmentation functions. The 
main tool in this respect is the availability of data at large z, for both light and heavy 
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nuclei, and for different hadrons. This combination explores the fragmentation process in 
a region which has not been investigated before. Obviously, our model has problems for 
z > 0.9 to explain the data. For such large z L, higher-twist processes in fragmentation, 
like the direct ejection of preformed pions in the nucleon, may be important. We think in 
other regions the data do not give a specific hint on the importance of higher-twist effects. 
Also it may be interesting to extend the predictions of other models of fragmentation to 
the full data set. 


6 Charge and flavour dependence 

The complete particle identification of the HERMES experiment allows us to disentangle 
the information for different hadron types. Recently, the HERMES collaboration reported 
measurement of R h M for pions, kaons, protons and antiprotons where the medium effects 
on 7r + and are equal, while significant differences are observed between K + and K~ 
and especially between p and p. The results for individual hadron species may reveal 
differences in the modification of q and q fragmentation functions as also suggested in 
Ref. [17], in the formation time of baryons and mesons [3,33], or in the size of the hadronic 
and prehadronic cross sections. 

In the following, we present the model results at the HERMES kinematics for charged 
pions and kaons, as the present version of this model does not apply to baryon production. 
The predictions are presented in Fig. 10, and compared with the HERMES data, both as 
a function of v and z. In each plot we show the contribution of the rescaling and the total 
contribution of rescaling and absorption for both charged states. One sees that the average 
charged meson spectra |(vr + + tt~) and \{K + + K~) are reduced by the rescaling of the 
singlet fragmentation functions, which are strongly affected by the anomalous dimension 
of gluon splitting. At large z, the observed ~10% effect can be calculated analytically. The 
difference of ( 7t~-tt + ) and (K~-K + ) spectra is governed by the non-singlet fragmentation 
function, which is weighted with the excess of neutron number (N=48) compared with 
the proton number (Z=36) in krypton. This explain the slight enhancement of n~ and 
K~ spectra at large z in the rescaling model shown in Fig. 10 (upper curves in each plot). 
Subsequent absorption with a larger nuclear cross sections for 7and K~ (see Fig. 7) 
corrects this difference. 

The data for 7r + and n~ production indeed show no difference in the multiplicity ratio, 
and are in nice agreement with the model predictions. We do not show calculations of 7r° 
mesons since they are very similar to the curve for charged pions. It is worth to point out 
that there is a nice agreement between model predictions and pion data also in the low z 
region. This Ending suggests that the observed difference in the spectra of charged hadrons 
(see Fig. 9) should be ascribed to proton production from target fragmentation which has 
not been considered theoretically. It is also possible to disentangle the contributions of 
the prehadronic and hadronic absorption by setting 0 y=O and <r* = a n in the calculation 
of R? M . The resulting ^-distribution of charged pions is shown as dotted curve in the 
upper left plot of Fig. 10. The neglect of hadronic absorption leads to a very small 
increase of R,% f at small z. This result shows that in the HERMES kinematic region 
mainly the prehadron absorption contributes because the hadron average formation length 
{h) = {^F) + f^ is larger than the size of the krypton nucleus. 
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Figure 10: Charge- and flavour-separated theoretical multiplicity ratios Rm(v) and Rm{z), com¬ 
pared with preliminary HERMES data on Krypton target [37] (the data for negatively charged 
hadrons have been slightly shifted to the right to improve the readability of the figure). The upper 
pair of curves includes rescaling without absorption for positive and negative particles, and the 
lower pair rescaling plus BC absorption for positive and negative particles. The dotted curve in the 
upper left plot show the result by setting <7^=0 and <T*=av in the calculation of R.J r . The error bars 
represent the statistical uncertainty only. Systematic uncertainty is about 4% [37]. 


For K + production there is also a general agreement between data and model compu¬ 
tations, while experimental results point to a stronger than calculated absorption of K~ 
mesons. We found that an increase of <r* of at least 50% would be necessary in order to 
reproduce the experimental data for K~. The discrepancy between theory and data may 
also point to a different formation mechanism for the negative kaons, as they do not con¬ 
tain any nucleon valence quarks, which dominate in the HERMES kinematics. This may 
imply a shorter formation length than predicted by Eq. (3.4), hence a larger absorption, 
since rank 1 hadrons would not participate in the K~ formation, see Eq. (B.16). 

The forthcoming HERMES data for different hadron types on light and heavy nuclei 
may help to further disentangle rescaling and absorption effects and to clarify details 
of the fragmentation process. In Fig. 11 we show the predictions of the model for Ne 
and Xe nuclei in the HERMES kinematic region for charged pions and kaons. To this 
respect, the absorption mechanism can be studied on heavy nuclei, while light nuclei are 
more sensitive to the rescaling effect alone. One characteristic feature of the rescaling 
corrections is the increase of the multiplicity ratio Rm > 1 for z < 0.2 which in the case 
of K + ,K~ fragmentation is visible in the calculation with rescaling alone as shown in 
Fig. 10. In addition, due to the smaller /F + -nucleus interaction cross section, the K + 
production on light nuclei may leave a chance to observe the pure rescaling effect. 
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Figure 11: Predictions of the rescaling plus absorption effects at HERMES for charge- and flavour- 
separated multiplicity ratios on Neon and Xenon targets. 

7 Discussion and Conclusions 

We have considered two effects to describe semi-inclusive hadron production in nuclei: 
gluon radiation and nuclear absorption. In the final theoretical result the two effects add 
and until now the experimental situation does not clearly identify the rescaling effect. 
Theoretically, the gluon radiation occurs through the whole hadronization process even 
outside of the nucleus and criticism may be raised whether the rescaling correction is 
inappropriate in this case. We do not believe so. Rescaling extends the gluon window 
in Q 2 toward the infrared and includes thereby those gluons which are later absorbed by 
target fragments related to the nucleus. The fast forward gluons which go into the current 
fragmentation region have a higher time-like virtuality and are taken into account in the 
process of absorption of the prehadron. Therefore we think that no special necessity arises 
to truncate the rescaling correction due to formation outside of the nucleus. Of course, 
by assuming only a constant rescaling factor over the whole nucleus, we do not consider 
the escape of quarks in surface nucleons directly into the vacuum. In this case a more 
sophisticated treatment may be needed. 

We have shown that rescaling models together with nuclear absorption are able to 
describe both HERMES and EMC data on the nuclear modification of hadron production 
in DIS process. While maximal deconfinement is ruled out by the data - it assumes a too 
large deconfinement - partial deconfinement is shown to have a possible effect on fragmen¬ 
tation functions. Nuclear absorption of the preformed hadron is shown to significantly 
affect the hadron production on heavy nuclei in the kinematic region of the HERMES 
experiments. The gluon-bremsstrahlung model of Ref. [16], that shares some similarities 
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with the “outside-inside” picture of the Lund model, is also in good agreement with data 
at z > 0.5. Further precise data at moderate and high z/’s and for light and heavy targets 
are needed to disentangle rescaling and formation length effects and to establish rescaling 
as an independent mechanism on its own. 

Different scenarios for the space-time development of hadronization may be usefully 
tested against HERMES data. In our theoretical description they enter through the prob¬ 
ability distributions V *, and eventually Vh if one separates between prehadron and hadron, 
and through the corresponding formation lengths (/*) = (Ip) and ( lh) = {If) + zv / ka- 
It would be interesting to measure ^-distributions with low- and high-// cuts in order to 
study how much the prehadronic absorption differs from hadronic absorption as the ex¬ 
perimental data on K~ may suggest. This would allow to change the characteristic length 
scale L = v /ka of the hadronization process. The high-// data sample would be sensitive 
only on prehadron absorption, because the hadron is formed almost entirely outside the 
nucleus. On the other hand, in the low-// data sample the information obtained on the 
prehadron absorption may be used to pin down the hadron absorption. Furthermore, the 
use of different cuts on z would allow to map out the //-dependence of the absorption pro¬ 
cess more neatly. The necessity to have very accurate data in a wide range of v and with 
different nuclei has also been stressed in the recommendation letter for a new European 
electron facility [39]. 
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Appendices 

A EMC and HERMES acceptance 

In Table 3 we give the list of the EMC and HERMES kinematic range and experimental 
cuts, used in the numerical computations of this paper. The HERMES kinematic accep¬ 
tance in z and v for pions and kaons is shown in Fig. 12. This figure explains the step 
structures of the LO computation of {v)(z) and (z)(u) shown in Fig. 8. 



EMC 

HERMES 


h 

h 

7r, K 

-E'beam 

GeV 

100 

200 

27.5 

27.5 

Qmin 

GeV 2 

2 

2 

1 

1 

W 2 - 

mm 

GeV 2 

4 

4 

4 

4 

I/max 


0.85 

0.85 

0.85 

0.85 

^min 


0.02 

0.02 

0.06 

0.02 

^max 


1 

1 

1 

1 

^rnin 


0.2 

0.2 

0.2 

0.2 

^max 


1 

1 

1 

1 

^min 

GeV 

10 

30 

7 

7 

^max 

GeV 

85 

170 

23.4 

23.4 

Zimin 

GeV 

3 

3 

1.4 

2.5 

-E'/imax 

GeV 

85 

170 

23.4 

15.0 


Table 3: Kinematic cuts of the EMC and HERMES experiments. 


pions & kaons 



Figure 12: The HERMES kinematic acceptance in z and v for pions and kaons. 
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B Lund formation lengths 


We review the probability distribution of the formation lengths of a hadron with fractional 
momentum z in the Lund model [12]: 


V h {y,z,L) 


Eg* Pn(y,z,L) 

Eg=i/ d yD n (,y;z,L) ’ 


(B.l) 


where D n (y; z, L) is the longitudinal phase space density for the n-th rank hadron with 
fractional momentum z to be formed at position y. For a prehadron with fractional 
momentum z, the probability distribution of the formation length is: 


'P,(y\z,L) = V h (y + zL\z,L) . 


(B.2) 


Consider Fig. 4. The n-th pair production point C n has the light cone coordinates 
Un ~ bi ± Vn■ The n-th rank hadron has light-cone momenta p± = E n ± p n , and the 
fractional momenta are: 


ZU Er 


Pn_ 

Po 


(B.3) 


The hrst constituent of the hadron h n , formed in C n , accelerates to the right under a 
constant string tension ka- At P n , it has a light-cone momentum p + = 2 ka (y n - i — Vn)- 
The second constituent of h n , created in C n _i, accelerates to the left, and at P n it has 

p + = 0. 


pi =ka (ri_i - yi) 


(B.4) 

(B.5) 


The hadron fractional momentum z n for n > 1 is determined by the coordinates of the 
pair creation points (recall ka/{2,v ) = 1/(2 L) and y ( | = 2L,z 0 = 1): 


vi~ i - vi 



(B.6) 


The probability that a quark fragments into n hadrons with fractional momenta z i, Z 2 , ... z n 
is related to the normalized string fragmentation function f(z) which gives the probability 
that a string breaks into two pieces carrying fractional momenta z and (1 — z) 


n n 

f(zi, z 2 ,..., Zn) n^=n ’ 

2—1 2—1 


with relative fractional momenta U; 


— Z n 


i ' n— 1 

1 _ Z2k =1 


From Eq. (B.6) for yf and using Eq. (B.7) we obtain: 


Vn = 2L ( 1 — ^2 Zk ) =2 L JJ(1 - u k ) 
k =1 / k =1 


(B.7) 
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The longitudinal phase space densities Di(y; z, L ) of the first rank hadron and D n (y; z, L ) 
of the n > 2-th rank hadron to be formed at y with fractional momentum z are: 

Di(y’, z, L) = J duif(ui)8(z — zi)S(y — L) (B.8) 

r n \ 

D n (y;z,L) = Y[du k f(u k )S(z-z n )6(y--[y+_ 1 -y-\) . (B.9) 

d k =1 

Since the n-th hadron production happens near the positive light-cone y~ ~ 0, we can 
use 

- n— 1 

5(z - z n )5(y - - [j/+_j - y~\ =S(z- u n y/L)S(y - L JJ(1 - u k )) 

k= 1 

to integrate D n over u n in Eq. (B.9). We obtain: 

D n (y, z, L ) = -f(—) Pn (y, L) , (B.10) 

y y 

where 

Pi(y;L) = 5{y — L) (B.ll) 

/ n n— 1 

Y[du k f(u k )S(y - L JJ(l-nj)) (B.12) 

k= 1 fc=l 

are the densities of formation points of the nth-rank hadron. The total density of hadron 
formation points is: 

OO OO 

ph(v, l ) = Y1 pn ^ h L ) = - L ) + pn ^ y ' L "> ■ ( B - 13 ) 

n =1 n =2 

Consider now the string fragmentation function of the Lund model [29]: 

/W=Ep«( 1 + C '«)( 1 -‘) C> . (B.14) 

<2 

where p q is a probability that a given production point involves pair production of a 
particular flavor q, and C q are parameters controlling the typical momentum fraction 
carried by final hadrons containing that flavor: 


C q 


C q = u,d, s 
D q q ^ u,d,s 


In the standard Lund model C ~ 0.3 and D q is flavour dependent. Note that the subscript 
q includes not only quark flavours, but also diquarks in the fragmentation of the string 
into baryons. Since p q •Cl for non light flavours q ^ u,d, s, we obtain from Eq. (B.12) 
and Eq. (B.13): 

/ T \ cv T \ IT c 

Ph(y, L)faS(y-L) T- . 

y 
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Finally, we obtain the probability distributions of the formation length of the prehadronic 
state, which for a flavour q = u,d,s has C q = C and for q ^ u,d, s has C q = D q : 


V2{v,z,l) 


1 + c q zL r y 1 9 

(1 + C)( 1 + Zj^) y~zL [(y + zL)( 1 - z)_ 

x j<% - (i - z)l\ + y^ C L d i( i - z ) L - y]| %] • 


(B.15) 


Note that the 5 -function represents the contribution from rank-1 hadrons. The average for¬ 
mation length for the prehadronic state is then computed as (l*)(z, L) = f dyyV*(y, z, L): 


(Q = 


1 + Cq 


1 + C+(C q -C)z 

1 + C (1 -z) 


x 


1 + 


2 + C 0 


y2-\~Cq 


2 -fi (2 + CL 2 + C q ; 3 + C g ; 


(1 — z) z L , (B.16) 


where 2^1 is the Gauss’ hypergeometric function [32], Subtracting 1 from the square 
brackets in Eq. (B.16) would give the average formation length for hadrons of rank n > 2. 
For a valence quark, q = u,d,s, the above expression reduces to Eq. (3.4), which is the 
one used in this paper. On the other hand, using Eq. (B.16) with C q =1 and 0=0 as 
considered in Ref. [34], the prehadron formation length is 


(Q 


ln(l /z 2 ) -1 + z 2 
1 — z 2 


zL 


5 


which is larger than in the standard Lund model as shown in Fig. 5. 


C Absorption formulae with two length scales 

When considering a different cross-section for the hadron and the prehadron interaction 
with the nucleons, a* 7 ^ cr^, equations (4.4) and (4.6) have to be modified to take into 
account the propagation of the prehadronic state. In the Lund model considered in this 
paper, it propagates through a fixed length equal to zv/ka before the formation of the 
hadron. 

In the Bialas-Gyulassy model Eq. (4.4) becomes 

ry'+zv/KA r °o 

W 0 (y',z,u) = 1 - 0 -* / dy*p A (b,y*) - a h dy*p A (b,y *) . (C.l) 

Jy' Jy'+zv/KA 

In the Bialas-Chmaj model, the two-stage fragmentation process is modeled as a chain 
decay with two average decay lengths: (/*) = (Ip) for the “decay” of the quark into the 
prehadronic state, and (A/) = zv/k for the “decay” of the prehadronic state into the 
observed hadron ((If) is given in Eq. (3.4)). Eq. (4.6) becomes then 

/»oo /»oo 

S A (b,y) = 1 - a* / dy'P*(y' — y)p A (b,y') — cr h / dy'P h (y'-y)p A (b,y') , (C.2) 

'by 
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where P*(y'—y) and Ph(y' — y) are, respectively, the probability of survival of a prehadron 
and a hadron at a distance y' — y from the virtual photon interaction point. When 
(/*) 7 ^ (A l) one has: 


P.W ~V)= (Q Al} m - e -(v'-»/(AQ) _ (C.3) 

W -y) = l - e-<»'-*'>«'•> - P,(y' - y) . (C.4) 

When (/*) = (A/) the prehadron survival probability is: 

P*(y' ~y) = i (y' - y )e~ {y '~ y)/{L} . (C.5) 
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